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The Microcanonical Fermionic Average method has been used so 
far in the context of lattice models with phase transitions at finite 
coupling. To test its applicability to Asymptotically Free theories, we 



D , have implemented it in QED2, i.e. the Schwinger Model. We exploit 



the possibility, intrinsic to this method, of studying the whole (3, m 
plane at negligible computer cost, to follow constant physics trajecto- 
ries and measure the m, — > limit of the chiral condensate. We recover 



H ' the continuum result within 3 decimal places. 



The Microcanonical Fermionic Average (M.F.A.) method for performing 
Lattice simulations with dynamical fermions fl]] is ideally suited for discussing 
the phase structure of theories with phase transitions at finite couplings, and 
it has been applied so far in this context 0, ^ . 

The conventional wisdom, however, requires that physically interesting 
theories are Asymptotically Free like QCD. It is then interesting to test the 
applicability of the M.F.A. method to a theory without phase transitions at 
finite coupling 0. In this paper we present an analysis of the Schwinger 
Model on the lattice. Strictly speaking, the Schwinger Model in the con- 
tinuum is not Asymptotically Free, since it is superrenormalizable and the 
Callan-Symanzik (3 function vanishes. However, in the lattice version, since 
the continuum coupling is dimensionful, the continuum theory is reached at 
infinite lattice coupling, much in the same way as four dimensional Asymp- 
totically Free theories like QCD. 

The continuum model is confining; it is exactly solvable at zero fermionic 
mass, so that we can compare the results of our simulations with exact ones. 
We have simulated the (unquenched) model in lattices ranging from 16^ to 
100^; we present here results for the average plaquette and for the chiral 
condensate, in the non symmetric {6 = 0) vacuum of the model. 

The evaluation of the chiral condensate has been made easier by the 
fact that, in the M.F.A. approach, the main computer cost resides in the 
evaluation of an effective fermionic action at fixed pure gauge energy by 
evaluating all the eigenvalues of the fermionic matrix at m = 0. It is then 
essentially possible, at no extra cost, to move in the plane /3,m to follow 
constant physics trajectories in approaching the correct continuum limit. 
This is easier in this model since here the Renormalization Group amounts 
to simple dimensional analysis. 

The M.F.A. method is fully described in [|1|. Starting from the partition 
function in terms of the total Action S = Sp + Sq, sum of the fermionic and 
pure gauge contributions, we define the density of states at fixed pure gauge 
Action (i.e. Euclidean Energy) 

N{E) = f DU6{Sg{U) - VE) (1) 

and an effective fermionic action through 

which is the microcanonical average of the fermionic determinant. 



In terms of the effective Action the partition function can thus be rewrit- 
ten as 

Z = f dEN{E)e-^^''-'^^ff^"''''f'''^ (3) 

Massless electrodynamics in 1 + 1 dimensions is confining, superrenormal- 
izable and exactly solvable. 
Its partition function is 

Z= f DA^DiPDijel'^'''^-^^^-'^''-'^^'^^^ (4) 

with the usual definitions of F^^, and D/ . The electric charge is 

The partition function (in the photonic sector) can be rewritten as |^ 

Z = /'Dv4^e-/''^'^[^^^'^^'"'+^^''^^l (5) 

i.e. as that of a theory of free massive vector bosons of mass M = -j^. In 
particular the Green's functions of purely bosonic operators are the same in 
both theories. This fact can be exploited for obtaining the average plaquette 
in the lattice (see later). 

As for the chiral properties of the model, the chiral current is anomalous. 
If the chiral limit is obtained from m 7^ 0, then the 6' = vacuum is selected. 
In this vacuum the chiral condensate is (with one flavour) 

-{^iP)c = 77^ = 0.15995 (6) 

while it diverges at zero flavour [i.e. the quenched limit) and is zero with two 
flavours. This is the value of the chiral condensate to be compared with the 
results of lattice simulations, where its chiral limit is obtained from m 7^ 

In the present simulation the pure gauge part is described in terms of 
non compact fields, while for the fermionic - gauge term we use the standard 
staggered formulation with n/ species. 

Since the photonic sector of the continuum theory is equivalent to a the- 
ory of a free, massive vector boson, the average plaquette of the Schwinger 
model can be compared with that of the vector boson, which can be exactly 
computed on a finite lattice: 

1 ^p 2-cospi-cosp4 .„x 

^ '"- 2V6:f .2dyJl- COS p^) + M^ ^^ 
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(P^ = I^^m) and for V^ -^ oo 

I r (Pp 2(l-cospi) + 2(1 -C0SP4) .„x 

The value M = ^ corresponds to the continuum Schwinger model, while 
the quenched value is 

E{M = 0) = ^ (9) 

Since Cc is dimensionful, (3 explicitely contains the lattice spacing: (3 = -^ 
so that the continuum limit of the theory is approached at (3 —>■ 00. The limit 
must be reached keeping fixed the dimensionless ratio ^^ = \/3m. This ratio 
defines constant physics trajectories. 

We have performed simulations in lattices up to 100^. We present here 
the results for the 64^ lattice, where we have the best statistics (for a total 
of 70 Cray-equivalent hours) 0. We will mainly discuss the 1— flavour case. 

As stated before, we compute all the eigenvalues of the fermionic matrix. 
This allows us to compute the Effective Action for all values of the mass, 
including m = 0. We have done so for 20 values of the energy, from 0.08 to 
1.3. 

One advantage of the MFA method is that the phase structure of the 
theory can be inspected directly from the fermionic effective action, whose 
derivatives must be discontinuous in order to generate a phase transition, at 
least for small nf, if the underlying pure gauge theory has no transition 0. In 
the case of QED2 the continuum theory is obtained as /3 — > cx) and one does 
not expects finite f3 transitions. The effective fermionic action numerically 
evaluated for the model does not show any sign of non analyticity and hence 
of phase transition p . 

The average plaquette is obtained as 

{E)l = ^^ ^ (10) 

and can be directly computed at m = 0. Since the underlying pure gauge 
theory is quadratic the density of states is known analytically 

N{E) = CaE^^-'i (11) 

so that the integrals in (10) are simple one- dimensional integrals. 

In Figure 1 we report the value of the average plaquette energy (dia- 
monds), multiplied by 2/3 to improve the visibility, compared with the exact 



result for a Massive Vector Model on the lattice. It is important to notice 
that the Schwinger Model is equivalent to a Vector Model in the continuum. 
On the lattice, there is no guarantee that the two models are related. From 
Figure 1 one can see that at small P, where presumably we are far from the 
continuum, there is disagreement between the numerical results and the an- 
alytical ones. However, already at /9 ~ 1 the agreement becomes excellent, 
showing that, at least for this operator, the continuum physics is reached 
fastly. The straight line is the quenched value 2f3{E) = 1, and one can see 
that, as m increases, the asymptotic value of (E) moves towards it. 
The chiral condensate 

cannot be directly computed at ttt, = 0, where it vanishes on the lattice, so it 
must be obtained as the limit m ^ 0. To reach the correct continuum value, 
this limit has to be taken simultaneously with the j3 -^ oo one, keeping the 
product \/]3m fixed. 

This can be easily done with this method, which does not require a sep- 
arate simulation of the fermionic contribution for each pair of parameters 
(Am). 

In Figure 2 we report the value of the chiral condensate for three values of 
the ratio ^^. For relatively large values of this ratio, scaling sets up already 
near /? ~ 1, but even for a very small value (in this lattice) as ^^ = 0.01, 
where finite spacing and volume effects appear in the small and large (3 
regions, there is a clear scaling window. We have repeated this procedure for 
12 values of ^^, and the values of the chiral condensate in the scaling window 
SO obtained have been reported in Figure 3. The behaviour of the condensate 
is very clear towards the continuum value, indicated in the figure as a circle. 
By fitting the points at small ^^ with polynomials we have always obtained 
consistent results for the intercept: 



(?/;^) = 0.160 ± 0.002 (13) 

in perfect agreement with the theoretical value. 

One can exploit the analogy between the fermionic system and a magnetic 
one to predict that the chiral condensate should behave as a power in ^^ in 
the scaling region, so allowing an unambiguous extrapolation at m = 0. 

From a formal point of view this result is also important, since it shows 
that, even with staggered fermions, where it cannot be proven rigorously, the 



usual introduction of the flavour number through powers of the ferniionic 
determinant is correct; in fact the numerical value for the chiral condensate, 
which exactly matches the continuum value, is obtained here by taking the 
square root of the determinant in the partition function. 

We have also analyzed the zero and two flavours cases |^. In the zero 
flavour limit there is really no scaling region, with the chiral condensate 
increasing at large (3, indicating that it diverges as expected 0. On the 
contrary in the two flavour case, the behaviour of the chiral condensate at 
finite mass indicates a vanishing value in the chiral limit, again in agreement 
with expectations ^. 

In conclusion, the results we find agree completely with the analytical 
expectations of the continuum theory. In this respect we believe that the 
MFA method can be applied to lattice models where the continuum limit is 
approached at infinite inverse coupling, like QCD. 

It is particularly interesting, in view of more ambitious applications, the 
ease with which constant physics trajectories can be followed in this ap- 
proach: in particular, since the mass dependence of the lattice Dirac operator 
has become trivial, it is possible to move in the (3, m parameter space at neg- 
ligible computer cost. It is useful to remember that also the Uf dependence 
is trivial 0. 

This potentiality has been fully exploited in the Schwinger Model, where 
Renormalization Group amounts to simple dimensional analysis and Con- 
stant Physics trajectories can be exactly defined through the whole param- 
eter space; as a consequence our numerical results for the chiral condensate 
are quite independent from the extrapolation to zero fermion mass and (as 
shown in Figure 3) are by far the best available in the literature. 

We believe that this potentiality can be used in more realistic theories 
like QCD. 

All the above simulations have been performed on various Transputer 
networks at L'Aquila University, Zaragoza University (RTN), the bulk on 
the Transputer Networks of the Theory Group of the Frascati National Lab- 
oratories of the INFN. 

This work has been partly supported through a CICYT (Spain) - INFN 
(Italy) collaboration. 
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Figure captions 



Figure 1 Average plaquette 64^, m = 0, nj = 1. 

Figure 2 Chiral condensate, 64^, ^ = 0.01, 0.04, 0.08 

Figure 3 Chiral condensate vs ^^, nj = 1, errors are smaller than symbols. 
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